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Abstract. An extension of the totaUy asymmetric exclusion process, which incorporates 
a dynamically extending lattice is explored. Although originally inspired as a model for 
filamentous fungal growth, here the dynamically extending exclusion process (DEEP) is 
studied in its own right, as a nontrivial addition to the class of nonequilibrium exclusion 
process models. Here we discuss various mean-field approximation schemes and elucidate 
the steady state behaviour of the model and its associated phase diagram. Of particular 
note is that the dynamics of the extending lattice leads to a new region in the phase diagram 
in which a shock discontinuity in the density travels forward with a velocity that is lower 
than the velocity of the tip of the lattice. Thus in this region the shock recedes from both 
boundaries. 
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1. Introduction 

When studying the physics of systems from the natural world, it is not surprising that one 
finds many problems which lie outside the realm of equilibrium physics. Rather, we are more 
likely to find systems which reach some sort of nonequilibrium steady state, characterised 
by non zero energy flow or other currents. For studying such systems, a vibrant collection of 
models has emerged. Many of the models have been studied extensively, not only for their 
applications but also for their intrinsically interesting behaviour [1]. 

The asymmetric simple exclusion process (ASEP) is a paradigmatic model for simple 
nonequilibrium driven diffusive systems [2]. It comprises a one dimensional lattice along 
which hard-core particles hop, with some bias in direction. Despite its simplicity, this 
model is capable of showing a range of interesting phenomena including boundary induced 
phase transitions, shock formation and spontaneous symmetry breaking. Over the years 
the apphcation of the ASEP has been wide and varied, but it particularly lends itself to 
the modelling of transport problems [3], for example the motion of molecular motors along 
microtubule filaments in biophysical systems IH El El [7j. Indeed the ASEP was originally 
introduced as a lattice model of ribosome motion along mRNA [8] . 

Recently, a particular case of the ASEP; the totally asymmetric exclusion process 
(TASEP), was generalised to incorporate a dynamically extending lattice [9l[10]. The new 
variation allows a particle reaching the end of a lattice to extend it by converting into a single 
new site. In this way, a connection is made between the microscopic constituents of the model 
and the dynamics of the system size. We also note that very recently an asymmetric exclusion 
process with Langmuir kinetics and a dynamic boundary has also been studied [11]. 

The model of [QKTO] was motivated by a problem in mycology: the process of filamentous 
fungal growth, where continuous growth of the filament tip is maintained by a supply of 
mass transported from behind the tip to the site of growth. In jSj it was argued that one 
could model the transport of vesicles by molecular motors along the network of microtubule 
fragments by an ensemble of continuous effective microtubules acting as dynamically 
extending asymmetric exclusion processes. Preliminary mean-field and simulation results 
showed that various growth regimes occurred, in particular regimes exhibiting a high density 
of motors near the tip. 

In this work, the extending TASEP of [9J is discussed in its own right, as a nontrivial 
addition to the existing class of exclusion process models. We discuss how steady states may 
arise in this growing system, and extend the well established TASEP phase diagram into 
a 3-dimensional parameter space. As with the TASEP, we may derive the phase diagram 
through a phenomenological approach which considers the dynamics of a shock that travels 
between the boundaries of the lattice. Through this approach, we arrive at a new subregion 
in the phase diagram in which a shock is moving away from both boundaries and thus 
persists in the system. We develop a series of mean-field approximations for our analysis. 
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Figure 1. Schematic of the open boundary TASEP with input rate a, hop rate p — 1 and 
extraction rate /3. 
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Figure 2. The phase diagram for the TASEP with schematic density profiles. Three distinct 
phases exist. For a or /3 < 0.5, two phases are possible; a high density phase for a > j3 and 
a low density phase for (3 > a. When a and (3 > 0.5, the system is in a phase of maximum 
current, no longer limited by the input or output rate. 

The simplest scheme, which for the ASEP actually predicts the correct phase diagram |12j. 
only predicts the phases qualitatively. However a refined mean-field approximation predicts 
phase boundaries which are rather close to those we determine by Monte Carlo simulations. 
A heuristic argument which considers the symmetry between particles entering one end of 
the lattice and 'holes' entering the opposite end provides a further prediction for the phase 
diagram which is very close indeed to that predicted by the refined mean-field theory. We 
find that this prediction is in fact the best fit to the simulation data. 

The paper is organised as follows: In section [2], we review the TASEP and introduce 
the dynamically extending model. In section [3] the model is formally defined and in section 
m a simple mean-field approximation is used to analytically approximate the steady states. 
In section |5] the mean-field approximation is extended to include correlations between sites 
1 and 2, resulting in a refined prediction of the model's phase boundaries. In section El we 
discuss the heuristic argument and the resulting phase diagram which appears to be in closest 
agreement with simulation. The Monte Carlo simulation results are discussed in section [71 
where phase diagrams for the model are also presented. Finally, we conclude in section [8l 



3 




a 



2 



3 



4 



N 



Figure 3. Schematic of the model with input rate a, hop rate p = 1, desorption rate /? and 
growth rate 7. 



2. A TASEP based model 

We begin by reviewing the well-known TASEP, upon which the dynamically extending 
exclusion process we study is based. It comprises particles that hop along a one-dimensional 
lattice in a single direction. No more than one particle may occupy each lattice site at any 
given time. Particles are injected at one end of the lattice and removed from the other end 
with rates a and (3, which are the only control parameters for the model. In the bulk of the 
lattice, the hops occur with rate 1. Figure [1] illustrates the TASEP schematically. 

The behaviour of this system is categorised in terms of its steady state phases. A steady 
state is reached when the density of each site, that is the time-averaged occupancy, no 
longer changes with time. Equivalently, when the system is in a steady state, the current of 
particles along the lattice is the same everywhere. The phases are defined by the system's 
macroscopic properties in the large system limit i.e. the particle density profile across the 
lattice and the particle current between any two neighbouring lattice sites. Phase transitions 
are induced by varying the boundary conditions, controlled by parameters a and p. The 
exact phase behaviour of the TASEP is well understood from the exact solution [I3l E] as 
well as mean-field [12] and other approaches p]. Three distinct steady state phases exist. 
For low input, a < 0.5, and a < (3, the system is in a low density phase, where the bulk 
of the lattice is at a low density equal to a and the current is equal to a (1 — a). When at 
a = (3, there is a discontinuous phase transition and the system enters a high density phase, 
now limited by a low output rate f3 < 0.5. The bulk of the lattice in this phase is at a high 
density equal to 1 — /? and the current is equal to — P). When a and j3 are both > 0.5, 
the system enters a phase of maximal current. The system is no longer controlled by the 
input and output rates, instead the bulk density is 1/2 and the current is 1/4. Furthermore, 
the density profile decays algebraically from the boundaries towards the bulk value 1/2. The 
results are summarised in the phase diagram of Figure El 

The model discussed in the following work has an additional third parameter, 7, which 
is the rate at which a particle reaching the ultimate lattice site may convert into a new 
lattice site, thus extending the length of the lattice by 1. We shall explore the resultant 
phase structure in the three dimensional parameter space of a, (3 and 7. 
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3. The DEEP 



The dynamicaUy extending exclusion process (DEEP) is defined by the rates at which the 
foUowing processes occur: particles in the bulk hop to the left with rate 1; particles enter the 
lattice from the right with rate a; at site 1 two processes may occur, particles detach from 
site 1 with rate (3 and particles may transform into a new lattice site with rate 7. These 
processes are illustrated schematically in Figure [31 Thus 7 is the parameter controlling the 
lattice growth, (3 allows particles to leave the end of the lattice without extending it and 
ratio 7//5 controls the efficiency with which the lattice extends. 

Indicating the presence of a particle by 1 and an empty lattice site by the stochastic 
dynamics is 

01 10 with rate 1 , (1) 
at site 1 1^0 with rate (5 (2) 

at site 1 1 ^ 00 with rate 7 , (3) 

at the rightmost site — 1 with rate a . (4) 

Note the direction of particle hopping is to the left (as opposed to most studies of the TASEP 
where the direction is to the right) for reasons which will soon become apparent. This choice 
is of course immaterial. 

Also note that under the simultaneous reversal of the direction of particle hopping and 
interchange of a and (3, the TASEP enjoys an exact particle-hole symmetry at the microscopic 
level, while the DEEP does not. This is due to 'growth' rule [3], which has no counterpart at 
the right boundary. 

3.1. Exact steady state equations 

In the DEEP model the system size is not conserved and continually increases, therefore 
care must be taken to define what is meant by a steady state. For example, the density at 
site i will generally only become stationary with respect to a particular frame of reference. 
Two natural frames to consider are the frame moving with the tip, in which case the tip is 
always at site 1 and the site label i measures the distance from tip, and the stationary frame 
where the right boundary is always at site 1 and the tip is at site where A^ increases as 
the lattice extends. A phase of the system is specified by stationary values of the densities 
and the frame of reference. Moreover we shall encounter one phase, the maximal current 
phase, that does not satisfy this criterion and instead is only quasi-stationary. 

To begin with, we choose to work in the reference frame of the growing lattice tip, so 
that the leftmost site is always labelled site 1. The input end of the lattice is site N{t) and 
process (HI) may be formulated as there being a reservoir of particles of density a at site 
N{t) + 1. As time t ^ 00 the system size N{t) —>■ 00, therefore we may write a boundary 
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condition at the input end of the lattice as 

a = hm (tn) . (5) 

We now construct the equations for the steady state dynamics. Note that when a growth 
event occurs, via process ([3]), all the previous lattice sites are relabelled i i + 1. The 
exact equations for correlations functions in the steady state are: 

d{n) 



dt 

d{T2) 

dt 
d{Ti) 



:0=([l-ri]r2)-(7 + /3)(ri), (6) 
^ = ([1 - t,]t,) - ([1 - ri]r2) - 7(rir2) , (7) 
= ([1 - r,]r,+i) - ([1 - Ti_,]Ti) + 7(rih_i - r,]) z > 3 . (8) 



dt 

Where is the occupancy at site i and can take on a value 1 or 0. The angle brackets denote 
an average in the steady state. Thus (rj) is the average occupancy in the steady state or 
particle density at site i. These equations may be understood in terms of a particle current 
along the lattice. The positive terms represent the current of particles entering the site from 
the right and the negative terms represent the current leaving to the left. The dynamics at 
sites 1 and 2 differ from those of the rest of the system due to the transition from particle to 
lattice site. In the steady state, the current 'in' balances the current 'out' everywhere. We 
thus obtain from ([MH]) the following exact expressions for the conserved steady state particle 
current J: 

J=ii + P){n), (9) 

J={r2[l-n]), (10) 
J= (r,+i[l-r,])-7(rir,) t>2. (11) 

The rate of lattice extension, or tip velocity, is given by: 

V = 7(ri) , (12) 

which leads us to a relation between the current and velocity: 

J = v{l + ^). (13) 

7 



Our approach to solve the system of equations flQ HTTl) is to take a mean field approximation in 
order to analyse the steady state equations and determine the phases and phase boundaries 
of the system. 

4. Simple mean-field theory 

In the simplest mean field approximation, we replace 

i^i^j) = PiPj V2,j(i^j) (14) 
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where = (rj) is the average steady state density at site i, and we ignore correlations 
between the density at different sites [12]. Equation (II ip then gives the recurrence relation: 

pi+i = -r-Z « > 2 . (15) 

-L Pi 

We find pi, p2 from (ED, ([lO]) and 

o - ^ - ^(^ + ^) (17) 

The recurrence f|T5|) together with the initial condition (|T71) determine all pi in terms of J 
(or t;). It remains to fix J self-consistently to determine the various phases. 

Low density Phase: The recurrence ( ITSl) has two fixed points 



P± = • (18) 

The lower of the fixed points is stable, so for p2 within the basin of attraction, the density 
iterates to this value and the density profile converges rapidly to a constant density in the 
bulk of the lattice. See Figure H] for a graphical representation of the iterative solutions to 
( IT5l) . To meet boundary condition ([5]), we define the lower fixed point: 

p- = a . (19) 

It follows from ( fTSi) that the upper fixed point is then: 

p+ = l-v-a . (20) 

The boundary condition ([5]) fixes v using f|T3l) and f|T5l) : 

7 

and the corresponding current is Jid = (1 + [3 /'~i)vid. 

The condition for this phase is that p2 < p+ so that p^ ^ p_ as z ^ oo. Using fl20l) and 
f l2T|) this condition reduces to 

This is the phase boundary for the low density phase. 

The density in the bulk of the low density phase is limited by a low input rate a. At 
the tip there is a decay to the boundary condition (|T9l) determined by rates j3 and 7. The 
low density region is split into two subregions distinguished by whether the tip density is 
above or below a. For p2 < a, the tip density is less than a and for a < p2 < 1 — — a, the 
tip density is greater than a. 
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Figure 4. A graphical representation of the solutions to the recurrence relation (|15p . 
Iterations beginning below the unstable upper fixed point, /Ci+, converge to the stable lower 
fixed point, 
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Figure 5. A shock is set up between regions in the lattice with a high density (next to 
the tip) and low density (far from the tip). The shock travels with a velocity fixed through 
mass conservation, which may be different to the velocity of the lattice tip. 



High density phase: Another possible solution to (fT5|) and ([5]) is that p2 = p+- In this case 
the recurrence (fTSll yields pi = p+ for i > 2. The phenomenological interpretation is that 
when 7, the rate of release of particles at the growing end, is no longer large enough to 
control the input rate; the particle density near the tip reaches a maximum value p+ that 
extends back through the lattice. The high density in the bulk p^ is then connected to the 
density at the input boundary a, with a discontinuous shock. Mass conservation implies a 
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speed for the shock, Vg, given by: 

Vs = 1 - p+ - a . (23) 

The shock solution is illustrated schematically in Figure [5l 

Defining p2 = p+, we obtain from f|T3|) . f|T71) and f|T8|) . the tip velocity in this phase: 

. 2(i_^ , (24) 

1-7 



and from this and (l20l) . we have an expression for the density at site 2 and throughout the 
rest of the high density region: 

p+ = l-27-/3. (25) 

From (122!) . one can see that this phase is entered when: 

7 = («-/5)/(l + «). (26) 

Due to the high density structure in the tip region we refer to the phase as the 'jammed' 
or high density phase. Note however that only on the phase boundary where Vg = v, is the 
shock and therefore the density profile stationary in the reference frame of the tip. In the 
rest of this phase Vs < v, and the shock moves away from the tip, resulting in an expanding 
high density region. Provided Vg > 0, the shock will also be moving away from the right 
boundary and the expanding high density region will never occupy the whole system. Using 
conditions (JTSl) and (|23i) . the condition Vs > becomes: 



1 + a 2 

In order to fulfil our criterion for a steady state (stationary densities in our frame of reference 

as t — s> oo) we should be in the reference frame of the shock. Then the density in front of 

the shock is given by 1 — 27 — /9 and the density behind is given by a. 

On the other hand if < 0, then the shock moves backwards through the system until 

it reaches the right boundary. In this region a high density is maintained throughout the 

length of the lattice except for a boundary region next to the right boundary. The system 

may be considered to be in a high density steady state with the reference frame being the 

stationary frame where the right hand boundary is fixed. The condition < becomes 

a — B , , 

7 < ^ • (28) 

Maximal current phase: Both the high and low density phases are bounded by the condition 
that the roots p± of ffT51) are real. Clearly, the roots ffTS]) are real provided (1 — vY > 4J, 
which implies using (fT3|) . that v < Vmax, where the maximum velocity satisfies 

2/5 + 37 - 2V/3^ + 37/3 + 27^ 



7 

and the corresponding maximum current is Jmax = {'^ + P / l)vn 



(29) 
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When J = Jmax, we have a single root of (fTS!) and from (fT8|) we deduce that the bulk 
density is 

1 —V 

'-'max ('iCW 

Pmax — 2 ■ v"-*^/ 

We refer to this phase as the maximal current phase; the flow is no longer controlled by 
the boundary condition a rather it has saturated at a maximal flow rate, Jmax- Note that 

Pmax < 1/2 and Jmax < 1/4. 

On a finite system of size N, in order to fix the boundary conditions, one takes 
J = Jmax + Oil/N"^) [121 [2]. The iterative solution of (ITSl) then yields a density profile 
decaying smoothly between both boundaries with an algebraic rather than an exponential 
decay of the density from the boundaries. Note that there is thus no natural reference frame 
in which we may define a steady state. Rather, the density at each site evolves as the system 
grows and we characterize this phase as quasi-stationary. We find the phase boundaries by 
considering the transitions to the maximal current phase from the low density phase, and 
from the high density phase: Jid = Jmax when 

and Jhd = Jmax when 

\ \ ■ (32) 
7-2 

3D Phase diagram: In order to visualise the phase structure of the DEEP, we may construct 
a 3-dimensional phase diagram from the results of the simple mean-field theory. Figure [6] is 
a schematic representation of the three phase regions, joined by the phase boundary surfaces 
described above. On the diagram, one can see that in the limit of no growth (7 = 0), the 
phase boundaries agree with those of the TASEP. 

4.1. The case /? = 

We may greatly simplify the equations of this model, by considering the case (3 = 0. This is 
the limit of no extraction from the tip and every particle reaching site 1 now contributes to 
the lattice growth. The current and velocity are therefore equal {J = v). We now have in 
the steady state: 

Pi = v/'j p2 = IV / {'J - v) (33) 

+ ^^2. (34) 

I- Pi 

The condition for the low density phase (12^ reduces to 7 > a/(l + a) and the tip 
velocity is given by 

v = J = ^^!r^. (35) 
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Figure 6. A schematic representation of the phase diagram in 3-D a, f3, 7 parameter space 
for the simple mean- field approximation. The three phases are shown (without subregions). 
The 7 = plane corresponds to the well known TASEP phase diagram. 



Then, the condition that p± be real becomes 

J = v <3- 2V2 , (36) 
which, with (1351) . imphes 

a<V2-l. (37) 
In the high density phase, where 7 < a/(l -|- a), the tip velocity (12^ is 

. = J = 2<1^. (38) 

For a /2 < '-f < a / {1 + a) a dynamic shock exists in the system. 

In the high density phase condition ( l36l) together with (l38l) imply 

7<l-^. (39) 

Therefore for the region bounded by a > \/2 — l, 7 > 1 — 1/ \/2, the system is in the maximal 
current phase, where the bulk density is pmax = — 1 and the current is Jmax = 3 — 2\/2. 

Monte-Carlo simulations of this system reveal that each of the predicted phases are 
indeed present. However, although the qualitative results of this simple mean-field theory 
appear to be correct, the predicted phase boundaries do not all agree with simulation (as we 
shall discuss in Section [7j). Thus in order to more accurately predict the phase boundaries, 
we develop a refined mean-field theory. 
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5. Refined mean-field theory 



In the approximation of the previous section, we ignored correlations between all pairs of 
sites in this system. In this section we introduce a refined mean field approximation which 
retains information about some correlations. 

We might expect that sites 1 and 2 are in fact strongly correlated because when site 
1 is vacated, so too is site 2. We can hope to improve our approximation by keeping the 
correlation {T1T2) intact in (jMlD and beginning the iteration of (fTS!) at site 3 instead of 2. 
Thus we approximate, for example, 

(riTar,) = (rirs)^ ^ > 3. (40) 

To begin the recurrence we must find a new expression for without factorising {T1T2), 
for example. Within the approximation, we have from (fTTi) . 

_ J + 7(rir2) 
1 - (ts) 

and from (fTOj) 

(rs) = J + (riTa) . 
Combining these two relations gives 

ir,r,) = PJ^ZA^ . (41) 

7 + P3 

Now we consider the steady state rates into and out of the configuration in which site 
1 is vacant and site 2 is occupied: 

"^^""'^^"^'^^ = = (rail - T2][l - n]) - {T,[l - n]) + P{t,T2) . (42) 
at 

We can identify the second term as the exact current from (fTOj) . Thus we obtain another 
exact expression for the steady state current: 

^= (r3[l-T-2][l-ri])+/5(rir2) . (43) 

Factorising beyond site 2, we may simplify and rearrange to give 

J - f3{TiT2) 

which along with (HTj) . gives us a quadratic equation for p3 

pI[1 - J - pA{i + m - J) - J] - {(^ + l)J = ^ ■ (45) 

This may be solved straightforwardly and used to deduce the three dimensional phase 
surfaces as in section HI 
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5.1. The case /? = 



We consider again the simple case of /5 = 0. Then (jHj) reduces to 

7 

The condition for a low density phase, ps < 1 — f — a, is now 

a < 7 (47) 

and the velocity as before is f = a(l — a)/(l + a). In the high density phase, 7 < a and the 
velocity is v = 7(1 — 7) /(I + 7). For a/ {2 — a) <7<aa dynamic shock exists between the 
high density tip region p = = (1 — 7)/(l + 7) and low density region p = a far from the 
tip. The maximum current phase is attained when a > — 1 and 7 > \/2 — 1. 



5.2. Further correlations 

We can repeat this mean-field improvement scheme by taking into account more correlations 
and correspondingly beginning the iteration of (ITSl) further from the tip. One might expect 
that as further correlations are taken into account, the predicted phase diagram will converge 
toward the exact solution. As we shall see however, this is not always the case. We do not 
present the details of the further refined mean-field theory here as the equations become 
rather cumbersome, but the results for the phase boundaries are shown in Figure [7] and 
discussed later, in Section [71 



6. Symmetry and heuristic argument 

An important feature which allows us to understand the phase behaviour of the TASEP 
(see Figure [2]) is the presence of a symmetry between particles hopping in one direction and 
'holes' hopping in the opposite direction. The high density particle phase can therefore be 
understood as a low density hole phase and vice versa. 

Due to the moving boundary condition and the transition ([3]), there is no exact 
microscopic symmetry between particles and holes in the DEEP dynamics. However, it 
is true that particles leave/holes enter at the tip whereas particles enter/holes leave at the 
right end. It is plausible then that there exists some symmetry between the high density 
phase (low density of holes) and the low density phase (high density of holes). Indeed, the 
phase boundaries between the high, low and maximal current phases predicted by the refined 
mean-field theory in the limit of /5 = display a symmetry under the interchange a ^ 7. 
But while the the velocities in the two phases, a(l — a.)/{l + a) and 7(1 — 7)/(l + 7), are 
also symmetric under the interchange a <-> 7 , the bulk densities in these phases, a and 
(1 — 7)/(l + 7), are not symmetric. 
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In the case of /3 7^ 0, the refined mean-field approximation no longer displays a particle 
hole symmetry. A heuristic argument however, does lead to a 3D phase transition between 
the high and low density phases which maintains a symmetry between input rate a and 
output rate /? + 7. 

If we ignore the density structure near the tip and assume that in the high density phase 
there exists a single high density value associated with the tip: p\ = pbuik, in the same way 
that in the low density phase, we have a single density a associated with the input end of 
the lattice, then from (llHllSp we have the following: 



I+P1 + 



This is solved to give the bulk density: 



1-7-/3 

Pi = Pbuik = — — , (49) 



1 + 7 

and associated high density velocity: 
7(1-7-/3) 



(50) 



1 + 7 

As before we appeal to the concept of a shock discontinuity in the system. The low 
density phase is reached when the shock velocity (1 — — a) is greater than the tip velocity 
f l50|l . which corresponds to a phase transition at: 

a = /3 + 7, (51) 

suggesting an exact correspondence between the balancing of input and output rates and the 
transition between high and low density phases. This boundary is the same as that derived 
by the refined mean-field theory in the limit f3 = and is in fact very close to that derived 
in the case /3 7^ (see Figure [TT] in the next Section). The remaining phase boundaries are 
derived as in section HI 



7. Phase Diagrams and Simulation 

As we have seen, the mean-field theories discussed above predict the same phases but with 
differing phase boundaries. The transition between the low density and maximal current 
phases ( 1311) is the same in each theory, however the transitions between other phases vary 
significantly. For comparison, the basic 2-D phase diagrams for the simple, refined and 
further refined mean-field theories in the limit of /3 = are plotted in Figure [3 Note in 
particular, a large difference in the boundary between subregions in the low density phase 
where the tip density decays from above or below a. The high to low density transition 
is also significantly different in the simple mean-field theory. In order to explore the exact 
phase structure of the DEEP model, and identify which of the theories best describes the 
exact behaviour of the model, we turn now to a numerical discussion. For simplicity, we 
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Figure 7. Comparison of three mean-field approximations in the case P — 0. 

work primarily in the (3 = limit. Note that in this limit the refined mean- field theory of 
Section [5] and the heuristic argument of Section [6] are in agreement. 

Monte-Carlo simulations of the DEEP were carried out across the parameter space. The 
simulated lattice comprises an array of binary values representing Tj. For each update, a site 
in the lattice is selected randomly and a transition is attempted with a probability defined 
by the rates (Pll]). One time-step consists of N such updates, so that on average each site is 
updated once per time-step. Thus as N increases, so does the number of updates per time 
step, in order to keep our time unit constant. Note that during one time-step the system 
may increase in length. In principle, this would affect our unit of time, however, such an 
increase, typically of one lattice unit, will be insignificant in the large N limit. 

Our aim is to allow the system to reach a steady state, and then we may calculate 
quantities such as site densities, particle fiux and rate of change in lattice length. A problem 
with this method is that as the system grows, N oo and the computational cost for each 
time-step increases. We may however minimise this problem by considering only a fixed 
length portion of the lattice in a particular frame of reference. For example, in the low 
density phase we expect that the density tends exponentially to a, so we may truncate the 
lattice at site M < N and simulate only the M sites nearest the tip. This method will only 
work for the steady state phases which are clearly defined in some reference frame. For the 
maximal current phase and shock region, we still consider the whole system. 

An example of each type of density profile is plotted in Figure El In the low density 
phase, the steady state densities are calculated by averaging the occupancy over many time- 
steps in the reference frame of the tip. On the other hand, in the high density phase, the 
steady states are defined in the reference frame of the input and the average occupancies 
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Figure 8. Density profiles from the different phases of the dynamicaUy extending exclusion 
process. From top left: Low density, shock profile, high density and maximal current phase. 
Note that since the high density profile is calculated in the reference frame of the input, the 
sites are labelled from right to left, site being where particles are injected. 



are calculated from this frame. For profiles in the shock region and the maximal current 
phase, the system is never in a truly stationary state, and to view the density profile, we 
must find the average occupancies at a fixed length and thus fixed time. The densities here 
are calculated by averaging over many arrays which have grown to a specified length. 

One property which may be used to locate the exact positions of the phase transitions 
is the lattice tip velocity, which is given by a simple function of the parameters and changes 
form across the phase boundaries (see equations ( 12T]|24|I29I) ). We find from the velocity of 
the simulated lattices, that the first refined mean-field theory appears to predict the correct 
form for the tip velocity and the correct phase transitions. Figure [9] shows the velocity of a 
simulated lattice plotted for different values of 7, against the refined mean- field expression 
for the velocity. The velocities closely follow the predicted curve f l35|) and as expected, a 
phase transition takes place at a = 7, where the velocity becomes independent of a. The 
maximum velocity is reached when a = a/2 — 1 and 7 > \/2 — 1. 

Further evidence in support of the refined mean-field theory can be found by exploring 
the density profiles of the low density region, where the three theories predict different 
transitions between profiles which decay at the tip from below and above a. In Figure [TOl 
we see that the transition appears to agree with that predicted by the first refined theory, 
while being at odds with both the further refined and the simple approximations. 

For completeness, we consider now the case of non zero j3. In Figure [TT] the simulated 
growth velocities are compared to the refined mean-field theory velocity functions for the 
high and low density phases for (3^0. The high density velocity obtained from the heuristic 
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Figure 9. The tip velocity plotted against a for 7 = 0.15. .0.45 and /? = 0. For 7 < \f2 — 1, 
a clear phase transition occurs when a ~ 7. This transition is from the low density phase 
where v = "'^l^^ to the high density phase, where the velocity is a independent. For clarity, 
dotted lines are used to indicate a = 7. When 7 = 0.45 > \/2 — 1, the maximum velocity 
occurs as predicted at a = — 1. 




Q. 



Ct site 

Figure 10. In the left-hand plot, the simple (solid line), refined (dashed line) and further 
refined (dotted line) phase boundaries are shown. Simulations were carried out in the low 
density phase, across the predicted boundary between profiles with a decay from above and 
below a (parameter values indicated with symbols). The resulting profiles, shown in the 
right-hand plot, display a transition that is in best agreement with the refined mean- field 
transition line. 
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Figure 11. Theoretical and simulation results for the tip velocity as a function of 7, given 
for fixed values of a and /3 > 0. The data points show simulation results, fitted to the 
heuristic and refined mean-field functions for the tip velocity in the high and low density 
phases. The low density velocity is the same for heuristic and mean-field approaches and is 
given by (|2ip . The high density velocity function is found in the refined mean- field theory 
by solving equation (|45p and inequality p3 < 1 — t; — a (explicit expression for the high 
density velocity is not presented here). The heuristic high density velocity is given by (|50p . 
The phase transition occurs at approximately vid = v^d 



argument is also shown and one may see that this coincides rather closely with the refined 
mean-field result. Due to the nature of Equations |45] and inequality ps < 1—v — a which must 
be solved simultaneously, in this case it is simpler to work with the velocity as a function 
of 7. Therefore, simulations were carried out with fixed values a = 0.4 and (3 = 0.2, in the 
range close to the phase transition 7 = 0.14. .0.28. We see a good agreement of the simulation 
data points with the high density and the low density velocity functions predicted by the 
refined mean field theory. Note however, that in the high density phase the simulation results 
appear after all to be in best agreement with the heuristic result for the velocity fl50l) . The 
phase transition occurs as expected where the high and low velocity functions meet. 

Finally, we present the detailed 2-D phase diagram of the first refined mean-field theory 
in the limit f3 = 0, shown in Figure [121 The diagram is split into the three distinct phases: 
high, low and maximal current. The low density phase is split into two subregions I and 
II, where the tip decays to the bulk density from above and below respectively. The shock 
region, where the shock recedes from both boundaries, is indicated within the high density 
phase. We believe this to be at least an accurate approximation to the exact phase structure 
as we discuss in the next section. 
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Figure 12. The phase diagram obtained from the refined mean-field approximation in the 
limit P = 0. Three phases are shown. The low density phase is split into subregions I and II, 
indicating regions where the density profile decays from above and below the bulk density 
respectively. The shock region is indicated within the high density phase. 

8. Conclusion 

In this work we liave studied a generalisation of the open boundary TASEP to a dynamically 
extending lattice and characterised a generalised phase diagram in the 7 parameter 
space. The model provides an example of a system where the dynamics of the microscopic 
constituents couple to and determine the dynamics of the containing structure. We conclude 
that steady states can exist in this system, where the lattice is continuously extending from 
one end while being supplied with particles from the other end. 

The well-known phases of the open boundary TASEP are reflected in the phase diagram 
of this model (Figure [T2|) . The low density phase has structure in the density profile near the 
tip and the bulk density is controlled by the low input rate at the opposite boundary. In the 
high density phase the bulk density is controlled by the tip dynamics and there is structure 
in the density profile near the opposite boundary. In between we have a shock region where a 
high density and a low density region are separated by a moving shock which in fact recedes 
from both boundaries. The shock region of the phase diagram extends the known phases 
and adds to the variety of shock formation phenomena which have been observed in the 
TASEP and related models[T5l[T6l[I7l[T8]. Finally, the maximal current phase is present but 
with bulk density < 1/2. However, in this model the maximal current phase is not strictly 
a steady state as the current and the density profile continue to evolve as the system grows. 

We have used several mean-field approximations and find that while the simplest 
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mean-field approach successfully predicts the qualitative behaviour of the phases, the phase 
boundaries are not accurately predicted. However, a refined mean-field approximation which 
takes into account the correlation between sites one and two predicts a phase diagram which 
appears to coincide very closely with simulation results and is very close to that deduced 
by the heuristic argument of Section [61 It would be of interest to know whether in fact 
the predictions, summarised for the P = case in Figure fT2l could be exact. Interestingly, 
taking into account a further correlation appeared to worsen the predicted phase diagram. 
This sort of negative result has been encountered in mean-field approximation schemes in 
the past [19]. For example it was found in the ZGB surface reaction model that a pair mean- 
field approximation successfully locates a discontinuous phase transition, while a four-site 
approximation seems to do worse |20j . 

The DEEP was initially formulated as a model derived from the TASEP of a biophysical 
growth process, namely fungal hyphal growth |9] and provides a simple framework which 
may find application in a broader context. The TASEP has been successfully developed to 
describe various biophysical transport problems by incorporating additional biological detail 
such as dynamic instabilities and heterogeneity [21], [221 EH] . It could prove fruitful to 
study the effects of such features in the present model. 
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